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A shear flow of particles in a laser-driven two-dimensional (2D) dusty plasma are observed in a 
further study of viscous heating and thermal conduction. Video imaging and particle tracking yields 
particle velocity data, which we convert into continuum data, presented as three spatial profiles: 
mean particle velocity (i.e., flow velocity), mean-square particle velocity, and mean-square fluctua- 
tions of particle velocity. These profiles and their derivatives allow a spatially-resolved determination 
of each term in the energy and momentum continuity equations, which we use for two purposes. 
First, by balancing these terms so that their sum (i.e., residual) is minimized while varying viscosity 
rj and thermal conductivity « as free parameters, we simultaneously obtain values for r\ and k in 
the same experiment. Second, by comparing the viscous heating and thermal conduction terms, we 
obtain a spatially-resolved characterization of the viscous heating. 

PACS numbers: 52.27.Lw, 52.27.Gr, 44.10. +i, 83.50.Ax 



I. INTRODUCTION 

Flows of most liquid substances are usually studied by 
modeling the liquid as a continuum, but there are some 
substances that allow the study of flows at the kinetic 
level, i.e., at the level of the individual constituent parti- 
cles. As examples, we can mention chute flows in granular 
materials [1] and capillary flows in colloids [2]. The solid 
particles in these soft materials are large enough that 
their motion can be tracked by video microscopy, allow- 
ing experimenters to record their positions and velocities. 
Like granular materials and colloids, dusty plasmas also 
allow direct observation of individual particle motion. 

Dusty plasma [3H2j is a four-component mixture con- 
sisting of micron-size particles of solid matter, neutral 
gas atoms, free electrons, and free positive ions. These 
particles of solid matter, which are referred to as "dust 
particles," gain a large negative charge Q, which is about 
— 10 4 elementary charges under typical laboratory condi- 
tions. The motion of the dust particles is dominated by 
electric forces, corresponding to the local electric field 
E = E con f + Ed, where E con f is due to confining poten- 
tials and Ed is due to Coulomb collisions with other dust 
particles. 

Due to their high charges, Coulomb collisions amongst 
dust particles have a dominant effect. The interaction 
force Q~Ed amongst dust particles is so strong that the 
dust particles do not move easily past one another, but 
instead self-organize and form a structure that is anal- 
ogous to that of atoms in a solid or liquid [8HT5]. In 
other words, the collection of dust particles is said to 
be a strongly-coupled plasma [TB] . In a strongly-coupled 
plasma, the pressure p is due mainly to intcrparticlc elec- 
tric forces, with only a small contribution from thermal 
motion [17] . 

Even when it has a solid structure, a collection of dust 
particles is still very soft, as characterized by a sound 
speed on the order of 1 cm/s [TSJ [T5]. As a result, a 



dusty plasma in a solid phase is very easily deformed by 
small disturbances, and it can be made to flow. Flows 
can be generated, for example, by applying shear using 
a laser beam that exerts a spatially-localized radiation 
force 19 27J . In such an experiment, the Reynolds num- 
ber is usually very low, typically R e rs 10, indicating that 
the flow is laminar [23j . 

This paper provides further analysis and details of the 
experiment that was reported in [25]. We now list some 
of the major points of these two papers, to indicate how 
they are related and how they differ. In this paper, we 
present: (1) a detailed treatment of the continuity equa- 
tions for both momentum and energy, (2) our method 
of simultaneously determining two transport coefficients 
(viscosity and thermal conductivity), (3) values of these 
two coefficients, and (4) spatially-resolved profiles of the 
terms of the energy equation, including the terms for vis- 
cous heating and thermal conduction, as determined by 
experimental measurements. In [25], we reported: (1) a 
discovery of peaks in a spatially-resolved measurement 
of kinetic temperature, (2) a demonstration that these 
peaks are due to viscous heating in a region of a highly 
sheared flow velocity, and (3) a quantification of the role 
of viscous heating, in competition with thermal conduc- 
tion, by reporting a dimensionless number of fluid me- 
chanics called the Brinkman number |29j which we found 
to have an unusually large value due to the extreme prop- 
erties of dusty plasma as compared to other substances. 
The values of viscosity and thermal conduction found in 
this paper are used as inputs for the calculations of di- 
mensionless numbers in |28] . The identification of viscous 
heating as the cause of the temperature peaks reported 
in |28j is supported by the spatially-resolved measure- 
ments reported here. 

In the experiment, the dust particles are electrically 
levitated and confined by the electric field in a sheath 
above a horizontal lower electrode in a radio-frequency 
(rf) plasma, forming a single layer of dust particles, 
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Fig. 1(a). The dust particles can move easily within their 
layer, but very little in the perpendicular direction, so 
that their motion is mainly two dimensional (2D). They 
interact with each other through a shielded Coulomb 
(Yukawa) potential, due to the screening provided by 
the surrounding free electrons and ions [3D]. As the dust 
particles move, they also experience frictional drag since 
individual dust particles in general move at velocities dif- 
ferent from those of the molecules of the neutral gas. This 
friction can be modeled as Epstein drag and char- 
acterized by the gas damping rate i/ gas . 

Using experimental measurements of their positions 
and velocities, the dust particles can of course be de- 
scribed in a particle paradigm. They can also be de- 
scribed by a continuum paradigm by averaging the par- 
ticle data on a spatial grid. In transport theory, mo- 
mentum and energy transport equations are expressed in 
a continuum paradigm, while transport coefficients such 
as viscosity and thermal conductivity are derived using 
the particle paradigm because these transport coefficients 
are due to collisions amongst individual particles. In our 
experiment, we average the data for particles, such as 
velocities, to obtain the spatial profiles for the contin- 
uum quantities, such as flow velocity. In the continuum 
paradigm, a substance obeys continuity equations that 
express the conservation of mass, momentum, and en- 
ergy. These continuity equations, which are also known 
as Navier-Stokes equations, characterize the transport of 
mass, momentum, and energy In a multi-phase or multi- 
component substance, these equations can be written 
separately for each component. The component of in- 
terest in this paper is dust. 

In Sec. II, we review the continuity equations for dusty 
plasmas. In Sees. Ill and IV, we provide details of our 
experiment and data analysis method. We designed our 
experiment to have significant flow velocity and signif- 
icant gradients in the flow velocity, i.e., velocity shear. 
In Sec. V, we simplify the continuity equations using the 
spatial symmetries and steady conditions of the experi- 
ment, and including the effects of external forces. We will 
use our experimental data as inputs in these simplified 
continuity equations in Sees. VI and VII. 



In this paper, the mass density p and the fluid velocity 
v describe the dust continuum. 

The momentum equation [3lJ E2] is 
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Here, p c , rj, and £ are the charge density, shear viscos- 
ity, and bulk viscosity, respectively; and rj/p is called the 
kinematic viscosity. In this paper, these parameters de- 
scribe the dust continuum. Equation |2]) describes the 
force per unit mass, i.e., acceleration, for the continuum. 
The confining field E con f and pressure p have been dis- 
cussed in Sec. I. The last term in Eq. ^ is due to the 
momentum contribution from forces such as gas friction, 
laser manipulation, ion drag, and any other forces that 
are external to the layer of dust particles, as discussed 
in Sec. IV. The other terms on the right-hand-side of 
Eq. ([2| correspond to viscous dissipation, which arises 
from Coulomb collisions amongst the charged dust par- 
ticles. The viscous term (?7/p)V 2 v was studied in a pre- 
vious dusty plasma experiment |23j . 

For our experiment, we can consider Eq. ^ for two 
conditions, with and without the application of the ex- 
ternal force f ox t. Without the external force, the fluid ve- 
locity v is zero so that Eq. ^ is reduced to p c E con f j p — 
S7p/p = 0, meaning that the confining electric force is 
in balance with the pressure (which mainly arises from 
interparticle electric forces). With the external force, the 
external confining force remains as it was without f ex t; 
moreover, the pressure p will be affected only weakly be- 
cause the density is unchanged (as we will show later) 
so that the interparticle electric forces that dominate 
the pressure will also be unchanged. Therefore, in us- 
ing Eq. ([2]), we will assume that the first two terms on 
the right hand side always cancel everywhere. 

The internal energy equation, as it is expressed com- 
monly in fluid dynamics |31l I32j , is 



T 



dt 



v Vs = $ 



(3) 



II. CONTINUITY EQUATIONS FOR DUSTY 
PLASMAS 



We now review the continuity equations for mass, mo- 
mentum, and energy for dusty plasmas. We will then 
discuss the significance of some of the terms for our ex- 
periments. 

The equation of mass continuity, i.e., conservation of 
mass, is 



dp 

~dt 



+ V • (pv) = 0. 



(1) 



Here, s is the entropy per unit mass, k is the thermal 
conducitivity, and T is the thermodynamic temperature 
of the dust continuum. The last term P GXt is due to the 
energy contribution from any external forces f ex t- 

We assume that the continuity equations, Eqs. (JlpJ), 
are valid for the dust particles separately from other com- 
ponents of dusty plasmas that occupy the same volume, 
such as neutral gas atoms. The coupling between the 
dust particles with other components is indicated in f oxt 
and P e xt, so that the momentum and energy of the dust 
particle motion is treated for the dust particles separately 
from other components. Other external forces, such as 
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those due to laser manipulation, are also indicated in f cxt 
and P xt- 

The first term on the right-hand-side of Eq. ^ is due 
to viscous heating. The viscous heating term $ depends 
on the square of the shear, i.e., the square of the gra- 
dient of flow velocity. A general expression for $ has 
many terms (cf. Eq. (3.4.5) of Ref. [SI] or Eq. (49.5) of 
Ref. |32j). but it can be simplified for our experiment by 
taking advantage of symmetries, as explained in Sec. IV. 

The second term on the right-hand-side of Eq. ([3| is 
due to thermal conduction. It arises from a tempera- 
ture gradient. Previous experiments with 2D dusty plas- 
mas include a study of the thermal conduction term in 
Eq. (§ [33]. 

In this paper, most of our attention will be devoted 
to the first two terms on the right-hand-side of Eq. Q. 
Using our experimental data, we will compare the mag- 
nitudes of these terms. In [28] . we demonstrated that 
viscous heating is measurable and significant, when eval- 
uated using only global measures like the Brinkman num- 
ber. Here we further evaluate viscous heating by charac- 
terizing it locally using spatially-resolved profiles for the 
terms in Eq. ([3|. We also develop a method to simulta- 
neously obtain values of two transport coefficients of rj 
and k. 

III. EXPERIMENT 

Here we provide a more detailed explanation of the ex- 
periment than in [28) . An argon plasma was generated in 
a vacuum chamber at 15.5 mTorr (or 2.07 Pa), powered 
by rf voltages at 13.56 MHz and 214 V peak-to-peak. 
We used the same chamber and electrodes as in [3"4] , 
The dust particles were 8.09 «m diameter melaminc- 
formaldehyde microspheres of mass irid = 4. 18 x 10" 13 kg. 
The dust particles settled in a single layer above the pow- 
ered lower electrode. The layer of dust particles had a 
circular boundary with a diameter of ~ 52 mm and con- 
tained ~ 10 4 dust particles. As individual dust particles 
moved about within their plane, they experienced a fric- 
tional damping [22j with a rate v gas = 2.7 s" 1 due to the 
surrounding argon gas. 

The particles were illuminated by a 488-nm argon laser 
beam that was dispersed to provide a thin horizontal 
sheet of light, Fig. 1(a). Using a cooled 14-bit digi- 
tal camera (PCO model 1600) viewing from above, we 
recorded the motion of individual dust particles. This 
top-view camera imaged a central portion of the dust 
layer, as sketched in Fig. 1(b). The movie is available 
for viewing in the Supplemental Material of [28] . The 
portion of the camera's field of view that we will analyze 
was 23.5 x 23.5 mm 2 , and it contained « 2500 parti- 
cles. We recorded > 5000 frames at a rate of 55 frame/s, 
with a resolution of 0.039 mm/pixel. Our choice of 55 
frame/s was sufficient for accurate measurement of var- 



ious dynamical quantities, including the kinetic temper- 
ature, although a slightly higher frame rate would have 
been optimal [35] . In addition to the top- view camera, 
we also operated a side-view camera to verify that there 
was no significant out-of-plane motion; this was due to a 
strong vertical confining electric field. Thus, we will an- 
alyze the particle motion data taking into consideration 
only the motion within a horizontal plane. 

At first, the dust particles self-organized in their plane 
to form a 2D crystalline lattice. The particle spac- 
ing, as characterized by a Wigner-Seitz radius [3B], was 
a = 0.26 mm, corresponding to a lattice constant b ~ 
0.50 mm, an areal number density no = 4.7 mm -2 , and 
a mass density p = nrrid — 1.97 x 10" 12 kg/mm 2 . Using 
the wave-spectra method for thermal motion of particles 
in the undisturbed lattice, we found the following param- 
eters for the dust layer: uj p d — 75 s _1 , Q/e — —9700, and 
cl/Xd = 0.5, where uj p d is the nominal 2D dusty plasma 
frequency [36] , Q is the particle charge, e is the elemen- 
tary charge, and Ad is the screening length of the Yukawa 
potential. 

We used laser manipulation [20Tf27] to generate sta- 
ble flows in the same 2D dusty plasma layer. A pair 
of continuous-wave 532-nm laser beams struck the layer 
at a 6° downward angle. In this manner, we applied 
radiation forces that pushed particles in the ±x direc- 
tions, as shown in Fig. 1. The power of each beam was 
2.28 W as measured inside the vacuum chamber. To gen- 
erate a wider flow than in [121 H3] ; the laser beams were 
rastcred in the both x and y directions in a Lissajous 
pattern as in [37], with frequencies of f x — 123.607 Hz 
and f y = 200 Hz. We chose these frequencies to be 
^> (jj p d/1iT to avoid exciting coherent waves. The rastered 
laser beams filled a rectangular region, which crossed the 
entire dust layer and beyond, as sketched in Fig. 1(b). 
This laser manipulation scheme resulted in a circulating 
flow pattern with three stable vortices, as sketched in 
Fig. 1(b). We designed the experiment so that in the 
analyzed region the flow is straight in the ±x directions, 
without any significant curvature. Curvature in the flow, 
which is necessary for the flow pattern to close on itself 
as sketched in Fig. 1(b), was limited in our experiment 
to the extremities of the dust layer, where it would not 
affect our observations. 



IV. PARTICLE AND CONTINUUM PARADIGMS 

We start our data analysis by analyzing data for indi- 
vidual particles, i.e., by working in the particle paradigm. 
Using image analysis software [35J, with a method opti- 
mized as in [39] to minimize measurement error, we iden- 
tify individual particles in each video image and calculate 
their x — y coordinates. We then track a dust particle be- 
tween two consecutive frames and calculate its velocity 
Vrf as the difference in its position divided by the time 
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FIG. 1: (Color online), (a) Side-view sketch of the appa- 
ratus, not to scale. A single layer of dust particles of charge 
Q and mass m,d are levitated against gravity by a vertical dc 
electric field. There is also a weaker radial dc electric field 
S con f which prevents the dust particles from escaping in the 
horizontal direction. Further details of the chamber are shown 
in [3l] . The two laser beams are rastered in the y direction so 
that they have a finite expanse, and they are offset in the y 
direction as shown in (b). (b) Top-view sketch of laser-driven 
flows in the 2D dusty plasma. In the region of interest, the 
flow is straight, with curvature of the flow limited to the ex- 
tremities of the dust layer. A video image of the dust particles 
within the region of interest is also shown. The region of in- 
terest is divided into 89 bins of width 0.26 mm so that particle 
data can be converted to continuum data. 



interval between frames [35] . Now having the position 
and velocity of all the particles in the analyzed region, 
we can study motion at the particle level. For example, 
in Sec. VI we will use data for the individual particles to 
calculate the rate of their energy dissipation due to their 
frictional drag with the neutral gas, P ex t- 

We next convert our data for individual particles 
to continuum data, i.e., we change from the particle 
paradigm to the continuum paradigm. This is done by 
averaging particle data within spatial regions of finite 
area, which we call bins. There are 89 bins, which are 
all long narrow rectangles aligned in the y direction, as 
shown in Fig. 1(b). Each bin contains « 30 dust parti- 



cles. We choose the shape of these bins to exploit the 
symmetry of the experiment, which in the analyzed re- 
gion has an ignorable coordinate, x. The width of each 
bin is the same as the Wigner-Seitz radius, a. To reduce 
the effect of particle discreteness as a particle crosses the 
boundaries between bins, we use the cloud-in-cell weight- 
ing method [TtJl HO] which has the effect of smoothing 
data so that a particle contributes its mass, momentum 
and energy mostly to the bin where it is currently located 
and to a lesser extent to the next nearest bin. The data 
are binned in this way regardless of their x positions, 
since x is treated as an ignorable coordinate. We also 
time-average these binned data, exploiting the steady 
conditions of the experiment. This procedure (binning, 
cloud-in-cell-weighting, and time averaging) yields our 
continuum quantities, such as the flow velocity v(j/). It 
also yields a kinetic temperature 
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fcs7k in = -m d \v d - v| : 



(4) 



which is calculated from the individual particle veloci- 
ties; this kinetic temperature is not necessarily identical 
to the thermodynamic temperature T. Here, fc^ is the 
Boltzmann constant. 

We assume that it is valid to use a continuum model 
when gradients are concentrated in a region as small as 
a few particle spacings. In fact, it has been shown ex- 
perimentally that the momentum equation |23) and the 
energy equation [33J remain useful in 2D dusty plasma 
experiments with gradients that are as strong as in our 
experiment. 

The notation we use in this paper distinguishes veloc- 
ities and other quantities according to whether they cor- 
respond to individual particles or continuum quantities. 
Parameters for individual dust particles are denoted by a 
subscript d, for example for the velocity of an individ- 
ual dust particle. Continuum quantities in a theoretical 
expression are indicated without any special notation, for 
example v for the hydrodynamic velocities in Eqs. 
and ([2]). Finally, continuum quantities that we compute 
with an input of experimental data, as described above, 
are indicated by a bar over the symbol, for example v. 



V. SIMPLIFICATION OF CONTINUITY 
EQUATIONS FOR DUST 

Here, we present our simplification of the continuity 
equations, Eqs. |IJ3]), to describe our 2D dust layer. 
These simplifications involve three approximations suit- 
able for the conditions in our 2D layer, and a treatment 
of two external forces, laser manipulation and gas fric- 
tion, that are responsible for f ext and P ex t- We describe 
these simplifications and the resulting continuity equa- 
tions, next. 
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A. Approximations to Simplify Continuity Equations 

The first of our three approximations is d/dt = 0. 
This approximation is suitable for the steady overall 
conditions of our experiment. Aside from the particle- 
level fluctuations that one desires to average away, when 
adopting a continuum model, the only time-dependent 
processes in the experiment were the rastering of the laser 
beam at > 100 Hz and the 13.56 MHz rf electric fields 
that powered the plasma. These frequencies are too high 
for the dust particles to respond, and the rastering of the 
beams is not a factor anyway because we will only use 
the continuum equations outside the laser beams. 

The second approximation is that d/dx is negligibly 
small. Due to the symmetry in our experiment design, 
as mentioned in Sec. Ill, x is treated as an ignorable 
coordinate, i.e., d/dx — 0. As a verification of this as- 
sumption, we observe that the ratio {dv^ / dx) / {dv^ / dy) 
is of order 10~ 2 , as a measure of the slightly imperfect 
symmetry of our experiment. Thus, for our experiment, 
we consider d/dy to be of zcroth order and d/dx to be 
two orders of magnitude smaller, when we approximate 
Eqs. ( [lP . 

The third approximation is that is negligibly small, 
based on our observation of the flow velocity of our 2D 
layer. Our results for the calculated flow velocity of the 
dust layer are shown in Fig. 2. From the velocity in 
Fig 2(a), the flow can be easily identified from the two 
peaks with broad edges. However, the flow velocity in 
the y direction is two orders of magnitude smaller than 
in the x direction, as shown in Fig. 2(b). 

Using the second and third approximations, and omit- 
ting terms that are small by at least two orders of magni- 
tude, we can easily see that we can approximate v • V = 
and V ■ v = 0. The latter indicates that the dust layer 
can be treated as an incompressible fluid in our experi- 
ment. Using these three approximations in Eq. 0, we 
find that Vp = 0. In other words, for our approxima- 
tions, the density p is uniform, which is confirmed by 
our experimental observation of Fig. 2(c). Because the 
density is so uniform, we can also assume that plasma 
parameters, such as the dust charge Q, are also spatially 
uniform within the analyzed region. 

In addition to these approximations, we also assume 
that r\ and k are valid transport coefficients for our sys- 
tem. There are theoretical reasons to question whether 
2D systems ever have valid transport coefficients, and 
this is typically tested theoretically using long-time tails 
in correlation functions [HJ 32] . It could be tested exper- 
imentally by repeating the determination of n and k for 
vastly different length scales for the gradients of velocity 
and temperature and verifying that they do not depend 
on the length scale. However, such a test is not practical 
for experiments like ours, which tend to have a limited 
range of diameters of dust layers that can be prepared. 
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FIG. 2: (Color online). Profiles of continuum parameters 
during laser manipulation, including (a-b) flow velocity, (c) 
areal number density, and (d) a measure \cf>s\ of local struc- 
tural order that would be a value of 1 for a perfect crystal [19] . 
Disorder, as indicated by a small value in (d), is found to 
be greatest where the shear is largest, not where the flow is 
fastest. 



B. External Forces 

Now we consider external forces that contribute mo- 
mentum and energy to the two equations, Eqs. ^ and 
(|]). For our experiment, we can mention six external 
forces acting on the dust layer: gas friction, laser manipu- 
lation, electric confining force, gravity, electric levitating 
force, and ion drag. Since we only study the 2D motion 
of dust particles within their plane, the last three forces 
are of no interest because they are in the perpendicular 
direction, and will not affect the horizontal motion of par- 
ticles that is of interest here. The electric confining force 
is balanced by the pressure inside the 2D dusty plasma 
lattice, /9 c E con f = Vp, as described in Sec. II. Thus, only 
two of the six forces need to be considered: gas friction 
and laser manipulation. 
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Gas friction is the main dissipation mechanism in our 
experiment. We can consider the effect of this friction 
first at the level of a single dust particle, and then at the 
level of a continuum. For the momentum equation, we 
note first that at the particle level, a single dust particle 
moving at a speed of experiences a drag acceleration 
of — ^gasVd- At the continuum level, the contribution of 
this drag to Eq. ^ is simply the average acceleration 
experienced by all dust particles in a given spatial re- 
gion, — t'gagV. For the energy equation, at the particle 
level the rate of energy dissipation for one dust particle 
is — 2i/ gas KE ( 2, which is the product of a drag force and ve- 
locity, where KE^ = is the kinetic energy of one 
dust particle. At the continuum level, averaging over all 
the dust particles in the given spatial region, the rate of 
energy loss per unit mass in Eq. Q3J) is —2i/ gas KE/ [15] . 



C. Simplified Continuity Equations 



transfer of momentum due to two mechanisms: viscosity 
arising from interparticle electric forces and frictional loss 
of momentum due to collisions with gas atoms. Equa- 



tion ( 10 1 describes the energy transferred from the vis- 



cous heating $ and thermal conduction (the second term) 
as being balanced by the energy dissipated due to fric- 
tion as expressed in the last term of Eq. ( |10| ). We will 
use Eqs. (9][T0) only in spatial regions where fi asC r = 
and Pi ascr = 0, i.e., outside the laser beam. In the next 
section, we will present our calculation of terms appear- 
ing m Eqs. ([6jl0| using the dust particles' position and 
velocity data from our experiment. The terms of inter- 
est in these equations are dvx/dy, d 2 v^/dy 2 , KE, and 
5 2 (v- v) 2 /dy 2 . 



VI. PROFILES OF QUANTITIES IN THE 
CONTINUITY EQUATIONS 



Using the three approximations listed above and taking 
into account gas friction and laser manipulation forces, 
the mass, momentum, and energy continuity equations 
become 



Vp = 



V d 2 v x 
p dy 2 



-+/* 



(5) 



(6) 



K Q2 

* + - tt^T - 2// gas KE/m d + P lascr = 0, (7) 
pdy 2 



where 



$ = 



V ( dv x 



p \dy 



(8) 



is the viscous heating term, after simplifications based on 
the assumptions of d/dx — 0, — 0, and V • v = 0. 

In this paper, we will restrict our analysis to a spa- 
tial region where the laser force and power are zero, i.e., 
fiascr = and Pi asC r = 0. In this region, the momentum 
and energy continuity equations will be further simpli- 
fied: 



d 2 v x 
dy 2 



4> 



p dy 



; T - 2^ gas KE/m d = 0. 



(9) 



(10) 



To simplify the problem, we will assume that r\ and n are 
independent of temperature, as discussed in [2"31 [55] , 

We can comment on the meaning of these two equa- 
tions. Equation ^ indicates a balance of the sideways 
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FIG. 3: (Color online). Profiles of the first (a) and second 
(b) derivatives of flow velocity. To make a comparison, the 
flow velocity profile is also provided in (c) . 

Our results for the first and second derivatives of the 
flow velocity are presented in Fig. 3. These results are 
calculated using the flow velocity profile in Fig. 2(a), 
which we reproduce in Fig. 3(c). The first and second 
derivatives in Fig. 3(a-b) will be used in Eqs. ^ and 
([9]), respectively. From Fig. 3(a), we can identify four 
points of maximum shear, i.e., maximum dv^/dy; these 
are at y = 4.1 mm, 8.3 mm, 13.6 mm, and 17.8 mm. 
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These points of maximum shear coincide with other fea- 
tures of interest: the minimum in the structural order, 
Fig. 2(d), and peaks in the mean-square velocity fluctu- 
ation profile, which we will present below. We find that 
disorder, as indicated by a small value in Fig. 2(d), is 
greatest where the shear is largest, not where the flow is 
fastest. 

Comparing panels (b) and (c) of Fig. 3, we find that 
the profiles for the flow velocity v x and its second deriva- 
tive are similar, in regions without laser manipulation, for 
example in the central region 8.3 mm < y < 13.6 mm. 
This similarity is expected from the momentum equa- 
tion, Eq. provided that the viscosity r\ is spatially 
uniform. We do not expect that r\ would be spatially uni- 
form since viscosity in general depends on temperature 
and the temperature is highly non-uniform, as we will 
show below. Nevertheless, we find that the two curves 
are nearly similar, with a small discrepancy that we will 
quantify below when we present data for the residual of 
Eq. j9). 
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FIG. 5: (Color online). Profiles of the mean-square particle 
velocity fluctuation shown separately for particle motion in 
the x and y directions, (a) and (b), respectively. These quan- 
tities combined correspond to the kinetic temperature T^n, 
Eq. Q. The profile of Tkm were reported in [28], where we 
discovered peaks in the kinetic temperature where the shear 
is largest. The second derivative of the averaged mean-square 
particle velocity fluctuation for the motion in the x and y di- 
rections (c) will be used to determine the second derivative of 
the thermodynamic temperature in Eq. (10 1. The thin curve 
in each panel is the flow velocity profile v^; its scale is shown 
in Fig. 3(c). 



FIG. 4: (Color online). Profiles of the mean squared particle 
velocity shown separately for particle motion in the x (a) and 
y (b) directions. These quantities correspond to the mean 
particle kinetic energy KE, i.e., (v^ + Vy)/2 — KE/m^. The 
thin curve in each panel is the flow velocity profile v^; its scale 
is shown in Fig. 3(c). 

Profiles of the mean squared particle velocity, corre- 
sponding to the kinetic energy in the energy equation, 



Eq. ( 10 ), are shown in Fig. 4. This kinetic energy includes 
energy associated with both the macroscopic flow and the 
fluctuations at the particle level. We will use these pro- 
files in determining KE in the next section, where we will 
find the residual of Eq. 



(10). 



In Fig. 5, we present our results for the mean-square ve- 
locity fluctuation, which corresponds to the kinetic tem- 



perature as in Eq. Q. We will use this kinetic temper- 
ature in place of the thermodynamic temperature T in 
the energy equation, Eq. (10). Unlike the kinetic energy 



KE, the kinetic temperature only includes the energy as- 
sociated with the fluctuations of particle velocity about 
the flow velocity v [44] . 

As reported in [28] , there are peaks in the kinetic tem- 
perature profile that coincide with the position of maxi- 
mum shear. These peaks can also be seen in Fig. 5(a-b). 
In [35], we attributed these peaks to viscous heating. As 
an intuitive explanation of viscous heating, consider that 
higher shear conditions lead to collisions of particles flow- 
ing at different speeds, causing scattering of momentum 
and energy that leads to higher random velocity fluctu- 
ations, and therefore higher kinetic temperature. In the 



next section, we provide further verification that the tem- 
perature peaks are due to viscous heating; we do this by 
confirming that three terms in the energy equation, in- 
cluding viscous heating, are in balance as indicated by 
their summing to zero. 



VII. RESULTS FOR RESIDUALS OF THE 
CONTINUITY EQUATIONS 
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We now examine the momentum and energy equations, 
Eqs. (|9| and (10), which are written so that the right- 
hand-side is zero. When we use these equations with an 
input of experimental data, however, the terms will not 
sum exactly to zero, but will instead sum to a finite resid- 
ual. We calculate these residuals, and we vary two free 
parameters, the viscosity r\ and the thermal conductivity 
k, to minimize the residuals. (Specifically, we minimize 
the square residual summed over all bins in the central 
region of 8.6 mm < y < 13.4 mm.) This minimization 
procedure yields the best estimation for the values for r\ 
and k, which will be our first chief result. We will then, 
as our second chief result, be able to make a spatially- 
resolved comparison of the magnitude of different terms 
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(b) 



in the energy equation, Eq. ( 10 ) 



Figure 6 shows the residual of the momentum equa- 
tion, Eq. ([9]). This result is shown for rj/p — 0.69 mm 2 /s, 
which is the best estimation of the kinematic viscosity. 
The data are shown as a spatial profile because we cal- 
culated Eq. ([£]) separately for each bin, i.e., each value 
of y. Two peaks in Fig. 6(a), located within the laser 
manipulation region, are due to the momentum contri- 
bution from the laser. We do not use Eq. ^ with these 
peaks because of a finite laser force fi aS cr there. Instead 
we will use the flatter region between these peaks, where 
fiascr = 0, as magnified in Fig. 6(b). The small resid- 
uals in this flatter region indicate that the momentum 
equation, Eq. ([9]), is able to accurately account for the 
momentum of our 2D dust layer, and that the minimiza- 
tion process in this region yields a value for the viscosity. 

Figure 7 shows the residual of the energy equation, 
Eq. (10 1, for n/(cp) = 8 mm 2 /s, which is the best es- 
timation of the thermal diffusivity. In this calculation, 
we used the value of r\j p = 0.69 mm 2 /s from the mo- 
mentum equation above, and we varied the value of n 
to minimize the residuals as described above. Two large 
negative peaks in Fig. 7(a) are due to the energy contri- 
bution from the laser manipulation Pi aS cr, which is not 
included in Eq. (10). The small values of residuals in 
the flatter region between these peaks, as magnified in 
Fig. 7(b), show that the energy equation, Eq. (10), ac- 



curately describe energy transport in our 2D dust layer, 
and that the minimization process in this region yields 
a value for the thermal diffusivity. By achieving a small 
value of the residual, we have verified that the three terms 
in the energy equation, Eq. (10), are in balance. Since 
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FIG. 6: (Color online), (a) A profile of the residual of the 
momentum equation, Eq. (|9J) , assuming a kinematic viscosity 
of 77/p = 0.69 mm 2 /s. For the central region, magnified in (b), 
the summation of the squared residual reaches its minimum 
when rj/p = 0.69 mm 2 /s; this minimization process is how 
we determine rj, which is one of our main results. The thin 
curve in (a) is the flow velocity profile v^; its scale is shown 
in Fig. 3(c). 



one of these terms is viscous heating and another is com- 
puted from the temperature profile which has peaks, the 
balance we observe here is consistent with the conclusion 
of [25] that the temperature peaks are due to viscous 
heating. 

As the first chief result of this paper, we obtain the 
kinematic viscosity value of v = rj/p = 0.69 mm 2 /s and 
thermal diffusivity value of K/(cp) — 8 mm 2 /s. These 
values are obtained simultaneously in a single experi- 
ment. A source of uncertainty in these values is system- 
atic error in u gas , for example due to particle size disper- 
sion or uncertainty in the Epstein drag coefficients |22) . 
This is so because our method actually yields results for 
rj/{pv s&s ) and n/(cpv gas ). Another source of uncertainty 
is our simplification that we neglect heating sources other 
than laser manipulation [43] ; in a test we determined that 
n/(cp) is in a range from 7.5 to 8.4 mm 2 /s, depending on 
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FIG. 7: (Color onlin e) . (a) A profile of the residual of the 



energy equation, Eq. (10 1, assuming a thermal diffusivity of 
X = ti/cp — 8 mm 2 /s. For the central region, magnified in (b), 
the summation of the squared residual reaches its minimum 
when n/cp — 8 mm 2 /s. This minimization process is how we 
determine k, which is another of our main results. The thin 
curve in (a) is the flow velocity profile v^; its scale is shown 
in Fig. 3(c). 



whether these small heating effects are accounted for. We 
note that our results of r\j p and k/ (cp) agree with the re- 
sults of previous experiments [231 [33] using the same size 
of dust particles and a similar value of a/Xo- 

Our experiment allows us to obtain spatially-resolved 
quantitative measurements of three heat transfer effects: 
viscous heating, thermal conduction, and dissipation due 
to gas friction (i.e., cooling). These three effects appear 
as the three terms on the left-hand-side of Eq. (10). As 



the second chief result of this paper, we plot these three 
terms, presented as spatial profiles, in Fig. 8. 

Examining the spatial profiles for these terms, in 
Fig. 8, we see the most prominent features are two large 
peaks for the gas dissipation term where the flow veloc- 
ity is fastest; in this region the energy dissipation due 
to gas friction reaches its maximum. Despite the promi- 
nence of these features, however, they are not what in- 
terest us here. Instead, we are more interested in the 
regions of high shear, near the edge of the laser manipu- 
lation. Recall that in these high shear regions, the flow 
velocity gradient dv^/dy is largest, and there are peaks 
in the profiles of the kinetic temperature and the second 
derivative of the flow velocity, Figs. 3 and 5. 

In Fig. 8, our spatially-resolved profiles reveal that vis- 
cous heating and thermal conduction terms are peaked in 
regions of high shear. The viscous heating term, Eq. Q, 
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FIG. 8: (Color online), (a) Profiles of three terms in Eq. ( 10 1, 
assuming the transport coefficients, rj and k, obtained above. 
For the central region without laser manipulation (b), thermal 
conduction is one order magnitude larger than the viscous 
heating. This is our second chief result: a spatially-resolved 
comparison of the different mechanisms for energy transfer in 
our 2D dust layer. The thin curve in (a) is the flow velocity 
profile its scale is shown in Fig. 3(c). 



is always positive, meaning that viscous dissipation is 
always a source of heat wherever it occurs. The ther- 
mal conduction term partial d 2 T/dy 2 , on the other hand, 
can be either positive or negative, indicating that heat is 
conducted toward or away from the point of interest, re- 
spectively. In locations where the shear is strongest, for 
example at y = 4.1 mm, the thermal conduction term 
is negative indicating that heat is conducted away from 
that point. 

Although viscous heating has great importance in all 
kinds of fluids, and it has been understood theoretically 
for a very long time [45], a spatially-resolved measure- 
ment of it is uncommon. In most physical systems vis- 
cous heating is usually hard to measure either because 
the temperature increase is overwhelmingly suppressed 
by rapid thermal conduction, as we discussed in |28j . or 



10 



because the thinness of the shear layer does not allow con- 
venient in-situ temperature measurements. Most experi- 
mental observations of temperature increases due to vis- 
cous heating are either external or global measurements, 
and not spatially-resolved measurements like those that 
we report here. Indeed, in our literature search, we 
found no previous spatially-resolved experimental mea- 
surements of the viscous heating term, not only for dusty 
plasma, but also for any other physical system. As we 
explained in [55], our ability to detect strong effects of 
viscous heating is due to the extreme properties of dusty 
plasma, as compared to other substances. Our ability to 
make spatially-resolved measurements is due to our use 
of video imaging of particle motion. 
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FIG. 9: (Color online). Profiles of (a) the ratio of the viscous 
heating and thermal conduction terms in Eq. ( 10 1, and (b) the 
absolute value of this ratio. The sign of d 2 T/dy 2 determines 
the sign of the ratio in (a) . A positive ratio indicates that heat 
is conducted toward the position of interest. The large values 
of the ratio in (b) at high shear regions indicate significant 
viscous heating at those locations. For comparison, in |28| 
we found a Brinkman number Br = 0.5 (indicated by the 
arrow) , which is a global measure of the flow that provides less 
detailed information than the spatially resolved ratio shown 
here. The thin curve in (b) is the flow velocity profile v^; its 
scale is shown in Fig. 3(c). 

To further analyze the second chief result of this pa- 
per, the spatial profiles of the viscous heating and ther- 
mal conduction terms in Fig. 8(a), we plot the ratio of 
these two terms in Fig. 9(a). This ratio has its largest 
positive and negative values in the regions of high shear. 
In Fig. 9(a), negative values of this ratio are observed 
to occur in the high shear regions, which indicates that 
heat is conducted away from these regions. This result is 
consistent with observation of kinetic temperature peaks 
here. To characterize the magnitude of these two terms, 



we plot in Fig. 9(b) the absolute value of this ratio. We 
can see that, within regions of high shear, this ratio can 
be as large as unity, or even larger. A typical value of 
this ratio in the shear region is of order 0.5 for our exper- 
iment. This matches the value of the Brinkman number, 
Br =0.5 that we found in [35J for the same experiment. 
The Brinkman number is a global measure of the viscous 
heating, in competition with thermal conduction. 



VIII. SUMMARY 

In summary, we reported further details of the laser- 
driven flow experiment in a dusty plasma that was first 
reported in |28j . We simplified the momentum and en- 
ergy continuity equations, exploiting the symmetry and 
steady conditions of the experiment. We developed a 
method to obtain transport coefficients by minimizing 
the residuals of continuity equations using the input of 
experimental data. As our first chief result, we use this 
method to simultaneously determine, from the same ex- 
periment, two transport coefficients: kinematic viscosity 
and thermal diffusivity, which are based on viscosity and 
thermal conductivity, respectively. As our second chief 
result, we obtained spatially-resolved measurements of 
various terms in the energy equation. We found that, in 
a laser-driven dusty plasma flow, viscous heating is sig- 
nificant in regions with high shear, which is consistent 
with the interpretation of [55J that the peaks in the tem- 
perature profile are due to viscous heating. 

This work was supported by NSF and NASA. 
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